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ABSTRACT 

Spectra of several R Coronae Borealis (RCB) stars at maximum light were ex- 
amined for the [Ci] 9850A and 8727A absorption lines. The 9850A line is variously 
blended with a Fe II and CN lines but positive identifications of the [C i] line are made 
for RCrB and SUTau. The 8727A line is detected in the spectrum of the five stars 
observed in this wavelength region. Carbon abundances are derived from the [C i] 
lines using the model atmospheres and atmospheric parameters used by Asplund et 
al. (2000). 

Although the observed strength of a Ci line is constant from cool to hot RCB 
stars, the strength is weaker than predicted by an amount equivalent to a factor of 
four reduction of a line's gj-value. Asplund et al. dubbed this 'the carbon problem' 
and discussed possible solutions. 

The [C i] 9850A line seen clearly in RCrB and SUTau confirms the magnitude 
of the carbon problem revealed by the C I lines. The [C i] 8727 A line measured in 
five stars shows an enhanced carbon problem. The g/-value required to fit the ob- 
served [C i] 8727 A line is a factor of f 5 less than the well-determined theoretical 
g/-value. We suggest that the carbon problem for all lines may be alleviated to some 
extent by a chromospheric-like temperature rise in these stars. The rise far exceeds 
that predicted by our non-LTE calculations, and requires a substantial deposition of 
mechanical energy. 

Key words: stars: abundances — stars: atmospheres — stars: variables — stars: R 
Coronae Borealis. 



1 INTRODUCTION 

The rare class of R Coronae Borealis (RCB) variables contin- 
ues to present fascinating puzzles. Two defining characteris- 
tics of these supergiants are the unpredictable fading of the 
star and the presence of helium, not hydrogen, as the most 
common constituent of the stellar atmosphere. The fadings 
are generally assumed to occur when a cloud of carbon soot 
forms on the line of sight to the star, but the formation 
process, and initial location and evolution of the cloud en- 
compass unanswered questions. How a normal star or a pair 
of stars evolves to possess a H-deficient atmosphere has chal- 
lenged theorists and observers alike. A leading idea is that an 
RCB star is the result of a merger of a He white dwarf with a 
C-O white dwarf. To develop observational tests of this and 



other ideas, one requires accurate estimates of the elemental 
abundances for a representative sample of the RCB stars. 

In attempting to supply this data, Asplund et al. (2000) 
uncovered yet another interesting puzzle which is discussed 
below. In RCB stars, as in normal stars, the strength of an 
absorption line is influenced by two principal factors: the 
ratio of the line to continuous opacity and the temperature 
gradient through the region of the atmosphere in which the 
line is formed. It has been suggested that photoionization 
of neutral carbon from high-lying levels appears to be the 
dominant source of continuous absorption across the visible 
spectrum (Searle 1961). A characteristic of RCB star spectra 
is a rich collection of absorption lines of neutral carbon from 
levels only slightly less excited than those contributing to the 
continuous opacity. Since the line and continuous opacity 
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originate from similar excited levels of the neutral carbon 
atoms, their opacity ratio is quite insensitive to fundamental 
parameters of the stars such as effective temperature, surface 
gravity, and the fractional carbon abundance (effectively, the 
C/He ratio once it exceeds 1% — see also Schonberner 1975; 
Cottrell & Lambert 1982). For weak lines, this insensitivity 
extends to differences in the turbulent velocities. 

This is observed in the RCB stars. A C I line retains 
its equivalent width even as 'metal' lines may vary consider- 
ably from star to star (see Figure 1 in Rao & Lambert 1996). 
However, the predicted equivalent widths of the C I lines ex- 
ceed the observed equivalent widths by a considerable factor. 
Expressed as the change in a line's (//-value needed to bring 
the predicted in tune with the observed equivalent width, 
the disagreement amounts to a factor of four. Asplund et al. 
dubbed this 'the carbon problem'. 

The principal purpose of this paper is to add a new clue 
to the solution to the carbon problem. That dimension is of- 
fered by the [C i] 8727 A and 9850 A lines and the estimate 
of the carbon abundance that they provide. We propose to 
investigate this aspect in more stars and to see whether a 
solution can be obtained for this puzzle. The present pa- 
per provides the results of these observations which, to our 
surprise, show an enhanced carbon problem. 



2 OBSERVATIONS 

High-resolution optical spectra of eight RCB stars (Table 
1) were obtained at the W. J. McDonald Observatory's 2.7- 
m telescope with the coude cross-dispersed echelle spectro- 
graph (Tull et al. 1995). In addition, 7 Cyg, a normal F-type 
supergiant was also observed (Table 1) at R = 120,000. The 
telluric absorption lines from the spectra of the programme 
stars were removed interactively using early-type rapidly ro- 
tating stars. We have used the Image Reduction and Analy- 
sis Facility (IRAF) software packages to reduce the spectra, 
and the task telluric within IRAF to remove the telluric 
absorption lines. 

The Th-Ar hollow cathode lamp fails to provide lines 
in the 9850A region for wavelength calibration. Plenty of 
Th-Ar lines are available shortward of 9670A to provide a 
wavelength solution to an accuracy of one tenth of a pixel. 
We apply this solution for the wavelengths longer to 9670A. 
To check the accuracy of the wavelength calibration in the 
9850A region, the measured wavelengths of the atmospheric 
H2O lines are compared with those measured from the solar 
spectrum (Swensson, Benedict, Delbouille & Roland 1970). 
We measured the wavelengths of 11 atmospheric H2O lines 
finding agreement to 0.003A ± 0.02A with those measured 
from the solar spectrum. 

Sample portions of R = 120,000 and, R = 60,000 spectra 
are shown in Figures 1, and 2. Note that, the spectra of the 
RCB stars are fully resolved at R = 60,000. All spectra are 
aligned to the rest wavelengths of N I lines which fall in the 
wavelength regions. Stellar lines were identified using the 
Revised Multiplet Table (Moore 1972), tables of spectra of 
H, C, N, and O (Moore 1993), Kurucz's list*, the Vienna 

* http://kurucz.harvard.edu 



Atomic Line Database^, the new Fel multiplet table (Nave 
et al. 1994), and a spectrum of 7 Cyg, a normal F-type 
supergiant, as a reference. 

All stars except V482 Cyg were observed at maximum 
light. Our spectrum of V482 Cyg was taken when the star 
was in a shallow minimum in July 1996, about 2 magnitudes 
below maximum light. It took the star about one year to 
return to maximum light. 



3 LINE DATA AND IDENTIFICATIONS 
3.1 The [Cl] lines 

The standard reference for the C 1 spectrum (Moore 1993) 
gives a predicted wavelength of 8727.126 A for the [Cl] 2p 2 
X D 2 - 2p 2 1 S transition and 9850.264 A for the [C 1] 2p 2 
3 P2 - 2p 2 X D2 transition. Wavelengths of the forbidden car- 
bon lines have not been measured directly from laboratory 
sources. 

The predicted wavelength of the 8727 A line is con- 
firmed by the wavelength of the absorption line in the solar 
spectrum (Allende Prieto, Lambert, & Asplund 2002) and 
the emission line in planetary nebulae (Liu et al. 1995). For 
the 8727 A line, we adopt loggf = -8.14 (Galavis, Mendoza 
& Zeippen 1997). The lower excitation potential is 1.264 eV. 
There is a blending Fe I line for which we adopt the log gf = 
-4.4 (Allende Prieto, Lambert, & Asplund 2002), but this 
line is not a significant contributor to the spectrum of an 
RCB star. 

For the 9850 A line, we adopt the predicted wavelength. 
There is supporting evidence from our spectra for this wave- 
length. The 2002 July spectra of RY Sgr show the 9850 A 
and 8727 A [Cl] lines in emission (see Figures 1 and 2). 
Adopting 8727.126 A as the rest wavelength for the latter 
emission and assuming the two lines have the same velocity, 
we find the astrophysical wavelength of the 9850 A emission 
is close to the predicted wavelength. Liu et al. (1995) report, 
however, a rest wavelength of 9850.36 ± O.OlA from detec- 
tions of the line in four planetary nebulae. The 0.1 A or 3km 
s _1 difference between the wavelength predicted from energy 
levels and that observed from planetary nebulae does not se- 
riously impact derivation of the carbon abundance from the 
line in the RCB star spectra. The transition probability for 
the 9850 A line is taken from Galavis, Mendoza & Zeippen 
(1997) who give A = 2.23 xl0~ 4 s" 1 or log gf = -10.79 
with an uncertainty of about 0.04 dex. The lower excitation 
potential is a mere 0.005 eV. 

Two additional lines complete multiplet IF. The 
stronger line with a p/-value about 3 times smaller than 
for the 9850 A line is at 9824.30 ± 0.05 A (Liu et al. 1995) 
and irretrievably blended in the spectra of cooler RCB stars 
(GUSgr, V482Cyg, SUTau, and RCrB). In the spectra of 
VZSgr, UVCas, and XX Cam, the 9824.30 A line is absent, 
but this is as expected from the S/N and the estimated 
strength of the observed 9850 A line. The third [Cl] line at 
9808.3 A is expected to be 4000 times weaker than the 9850 
A line, and undetectable. (The multiplet 2F 2p 2 3 Pi - 2p 2 
1 So transition at 4621.57 A with loggf = —11.12 is almost 

t http://www.astro.univic.ac.at/~vald 



Table 1. Tabic of observations; the stars are ordered by their effective 
temperatures decreasing from top to bottom 



Star 


Date 


R = 
A/AA 


S/N 
at 8727 A 


S/N 
at 9850 A 


Vrad 

km s _1 


RCB stars: 












XX Cam 


9 Oct f997 


60,000 




185 


11 


XX Cam 


26 Jan f998 


60,000 




120 


13 


XX Cam 


17 Nov 2002 


60,000 




186 


11 


XX Cam 


17 Nov 2002 


60,000 


278 






RY Sgr 


22 June 1997 


60,000 




170 


-30 


RY Sgr 


31 July 2002 


120,000 


300 


200 


-12 


UV Cas 


24 July 1996 


60,000 




140 


-32 


UV Cas 


15 Nov 2002 


60,000 




120 


-32 


UV Cas 


15 Nov 2002 


60,000 


178 






VZ Sgr 


18 Aug 1999 


60,000 




80 


245 


R CrB 


17 June 1995 


60,000 


340 


180 


20 


R CrB 


7 Aug 1995 


60,000 




110 


22 


R CrB 


20 June 1997 


60,000 




185 


21 


R CrB 


31 July 2002 


120,000 


400 


200 


19 


SU Tau 


15 Nov 2002 


60,000 




100 


43 


SU Tau 


15 Nov 2002 


60,000 


180 






V482 Cyg 


24 July 1996 


60,000 




90 


-40 


GU Sgr 


22 June 1997 


60,000 




95 


-45 


Standard star: 












7 Cyg 


31 July 2002 


120,000 


520 


230 





certainly blended as the blue spectral region of an RCB star 
is rich in strong lines.) 

3.2 Overview of the Spectra 

The region around 9850 A is shown in Figure 1 where the 
stars are ordered by increasing temperature from bottom to 
top. The spectrum of the normal yellow supergiant 7 Cyg is 
shown at the bottom. In the spectra of the hottest stars, the 
strongest lines are from C I and N I transitions with weaker 
lines of Si I, Fel, and Fell. The [Cl] line falls in the red 
wing of a Fe 11 line. Spectra of the coolest stars GU Sgr and 
V482 Cyg show additional lines which we identify as high 
rotational lines of the 1-0 band of the CN molecule's Red 
System (Davis & Phillips 1963). These CN lines are promi- 
nent in the spectra of GU Sgr and V482 Cyg and traceable in 
SUTau, RCrB and RYSgr. Two CN lines bracket the [Cl] 
line. 

Inspection of Figure 1 shows several additional features 
of interest. As anticipated, the strength of a given C I line is 
unchanged along the temperature sequence from XX Cam at 
the high temperature end to GU Sgr at the low temperature 
end. Also anticipated is the very large width of all absorption 
lines in the spectra of RCB stars. Compare the widths with 
those of lines in 7 Cyg, which itself is commonly referred to 
as showing broad lines. A surprise is the appearance of emis- 
sion in the RY Sgr spectrum near the wavelength of the [C 1] 
line. The illustrated spectrum is from 31 July 2002 in which 
the 8727 A line is also in emission. A spectrum from 22 June 
1997 also shows emission but to the red of the wavelength 
of the [C 1] line. The possibility of emission contaminating 
the absorption line is an unfortunate complication. 

The [C 1] 8727 A line (Figure 2) is slightly blended with 
a Si 1 line. Emission is very pronounced in the spectrum of 
RY Sgr. The central line depths of the [C 1] line run from 
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Figure 1. The spectra of the RCB stars in the region of [Cl] 
9850. 264A line, which is indicated by a vertical solid line, are 
plotted with their effective temperatures increasing from bottom 
to top. The spectrum of 7 Cyg is also plotted. The positions of a 
number of other atomic and some CN lines are also marked and 
the vertical dashed lines indicate the Fe II features. The relative 
intensity scale for GU Sgr and XX Cam is different from the other 
RCB stars. 

about 16 % for SUTau to 36 % for RYSgr with no obvious 
systematic trend with the effective temperature. 
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Figure 2. The spectra of RCB stars in the region of [Cl] 
8727. 126A line, which is indicated by a vertical solid line, arc 
plotted with their effective temperatures increasing from bottom 
to top. The spectrum of 7 Cyg is also plotted. The positions of 
the key lines are also marked. 

4 SPECTRUM SYNTHESIS - PROCEDURE 

The fact that the [C 1] lines are blended with other lines re- 
quires that the carbon abundance be extracted by spectrum 
synthesis. In addition to synthesizing the [Cl] line and its 
blends, we use lines of Cl, NT, Sil, Fel and Fen to obtain 
abundances and compare with values given by Asplund et 
al. (2000). 

Synthetic spectra are generated using the H-deficient 
model atmospheres computed by Asplund et al. (1997) and 
the Uppsala spectrum synthesis code BSYNRUN. Models 
corresponding to a C/He ratio of 1 % by number were cho- 
sen. This composition is equivalent to a carbon abundance 
loge(C) = 9.54 on the scale E^e; = 12.15, where /i; is the 
mean atomic weight of element i. Stellar parameters derived 
by Asplund et al. (2000) were adopted as our initial values. 
Synthetic spectra were convolved with a Gaussian profile to 
account for the combined broadening of stellar lines by the 
atmospheric macroturbulent velocity, rotational broadening, 
and the instrumental profile. 

4.1 The 9850 A Region 

The 9850 A [C 1] line may be blended with Fc n and CN 
lines. The Fell line at 9849.74 A has a lower excitation po- 
tential of 6.729 eV. To obtain the <?/-value of the line, we use 
the spectrum of 7 Cyg where the Fen line is present. With 
a model atmosphere constructed for the stellar parameters 
and iron abundance given by Luck & Lambert (1981), we 
find log gf = -2.60. The basic data for the CN 1-0 lines 
were taken from Kurucz's line list and checked against the 
SCAN tape (J0rgensen & Larsson 1990) and a recipe given 
by Bakker & Lambert (1998). A dissociation energy of 7.65 
eV was adopted for the molecule (Bauschlicher, Langhoff & 



Taylor 1988). Unblended CN lines away from the 9850 A 
feature were fitted by adjusting the N abundance, and this 
abundance used in computing the CN lines blended with the 
[C 1] line. An estimate of the N abundance was obtained not 
only from the CN lines but also from N I lines near the [C 1] 
lines. The gf- values for the N I lines were taken from Wiese, 
Fuhr & Deters (1996) compilation. 

The C 1 lines at 9852.27 A, a blend of two lines, and 
9859.15 A (Figure 1) have also been used for estimating the 
C abundance. Their <?/-values are determined by inversion 
of their equivalent widths in the 7 Cyg spectrum using the 
C abundance derived from a set of weak C I lines at shorter 
wavelengths. 

4.2 The 8727 A Region 

The red wing of the [C 1] line is blended with a Si I line 
(Figure 2). The gf- value of the latter line is derived from 
its equivalent width in the solar spectrum. The <?/-value of 
the Nl line to the red of the Sil line is taken from Wiese, 
Fuhr & Deters (1996). Spectrum synthesis of this region was 
confined to a 6 A window around the [C 1] line. 



5 SPECTRUM SYNTHESIS - RESULTS 

The [C 1] lines were sought in order to shed new light on 
the 'carbon problem'. Asplund et al. (2000) suggested that 
the 9850 A [C 1] line was not subject to the carbon problem. 
Their analysis was flawed in that the contributions of the 
blending Fe n and CN line were not considered. 

In addition to fitting the [C 1] lines, we analysed other 
lines in order to check the abundances given by Asplund et 
al. (2000). To within the uncertainties of the analyses and 
with the exception of V482 Cyg (see below), we confirm the 
result for carbon using C I lines in the 9850 A region, the N 
abundance using Nl lines in the 9850 A and 8727 A regions, 
the Si abundance from Sil lines in the 9850 A and 8727 A 
regions, and the Fe abundance from Fel and Fell lines in 
the 9850 A region. The CN 1-0 lines are also used to check 
the N abundance. In these atmospheres, molecules are trace 
species whose formation does not affect the partial pressures 
of the constituent atoms. Then, consideration of the ratio of 
line to continuous opacity shows that the strength of a CN 
line should be dependent on the N abundance but indepen- 
dent of the C abundance. 

Following remarks about 7 Cyg, we discuss the RCB 
stars in order of decreasing effective temperature. It should 
be noted that the C abundances given here for the RCB stars 
are not self-consistent - the model atmospheres and their 
continuous absorption were calculated assuming a logarith- 
mic C abundance of 9.54. If the same abundance is adopted 
for the calculations of the carbon lines - permitted and for- 
bidden - the predicted and observed equivalent widths are 
in disagreement. Thus, the C abundance derived here for the 
RCB stars is basically a measure of the inconsistency, which 
we term the carbon problem. 

5-1 7 Cyg 

Our analysis of permitted and forbidden C I lines uses a 
MARCS model atmosphere (Gustafsson et al. 1975) for 
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Figure 3. Observed (solid line) and synthetic [Cl] 9850. 264A line 
profiles of RCB stars. Synthetic profiles, including the blends, arc 
shown for each star for three different abundances. The maximum 
abundance of log e(C) = 9.5, the model input abundance is shown 
by the long dashes and in no case does this fit the observed fea- 
ture. The minimum abundance of loge(C) = 7.5 shows in effect 
the predicted spectrum for no contribution from the [C i] line (see 
dot-dash line). The best-fitting synthetic spectrum is shown by 
the line made of short dashes and the corresponding abundance 
is indicated in each panel. 



the atmospheric parameters derived by Luck & Lambert 
(1981). The 8727 A [Cl] line (Figure 2) yields the abun- 
dance log e(C) = 7.88. The 9850 A [C i] line is not detectable 
and the upper limit to its equivalent width corresponds to 
loge(C) < 8.18. A selection of 14 C I lines from 5380 A to 
8873 A with gf- values from Wiese, Fuhr & Deters (1996) 
give the abundance log e(C) = 7.94 ± 0.14. These lines range 
in equivalent width from 6 mA to 82 mA. For 7 Cyg, it is 
seen that the C I lines and the 8727 A forbidden line return 
the same abundance within the uncertainties. This result is 
in sharp contrast to the results from the RCB stars. 

5.2 XX Cam 

The 9850 A feature is dominated by the Fe n line such that 
we can set only an upper limit to the [C 1] contribution and 
the associated C abundance of loge(C) < 8.8 (Figure 3). 
The 8727 A line (Figure 4) is fit with an abundance log e(C) 
= 8.4. The limit from the 9850 A line to the C abundance 
is consistent with that from the 8727 A line. The two C I 
lines near 9850 A give an abundance consistent with the C 
abundance of log e(C) = 9.0 ± 0.4 from Asplund et al. who 
used a collection of C 1 lines and Opacity Project gf- values 
(Lou & Pradhan 1989; Hibbert et al. 1993; Seaton et al. 
1994). 

Asplund et al. put the uncertainty in T e g at about ± 250 
K which corresponds to corrections to abundances derived 
from C 1 and [C 1] lines of 0.05, and 0.2 dex, respectively. The 
uncertainty in logp ~ ± 0.5 provides only minor corrections 



Figure 4. Observed [Cl] 8727. 126A line profiles (solid lines) of 
RCB stars. Synthetic spectra are shown for four C abundances, 
as shown on the figure. 



to the abundances. An additional potential source of uncer- 
tainty is the adopted microturbulent velocity of £ = 9.0 km 
s _1 (Asplund et al. 2000). A change of £ by ± 2 km s _1 
changes the C abundance by 0.2 dex from the C I lines and 
by only 0.05 dex from the [C 1] 8727 A line. These estimates 
of uncertainty do not allow the 0.6 dex difference in the C I 
and [C 1] abundances. 

The carbon problem for the [C 1] line is even more severe 
as seen for 8727 A [C 1] line which gives an abundance 0.6 less 
than the C I lines and 1.1 dex less than the input abundance. 



5.3 RY Sgr 

The 31 July 2002 spectrum shows emission in the blue wing 
of the 8727 A forbidden line. Peak emission occurs at the 



velocity of —29 km 



This is blue-shifted by about 17 



km s _1 from the photospheric velocity of —12 km s _1 . The 
emission line has a width (FWHM) of about 19 km s~ . A 
spectrum taken the same night shows emission at 9850 A 
(Figure 1). Emission at [Cl] may be a common occurrence 
for this star. Our 1997 June 22 spectrum shows 9850 A in 
emission at the radial velocity of —10.5 km s _1 equivalent 
to a red-shift of 19.5 km s _1 relative to the photosphere. 
This spectrum does not include the region around the 8727 
A line. 

RY Sgr exhibits a pulsation with an amplitude of about 
35 km s" 1 (Lawson, Cottrell & Clark 1991). The mean or 
systemic velocity was given as —21 km s _1 by Lawson, Cot- 
trell & Clark (1991). Relative to this velocity the emission 
at 8727 A in the 2002 spectrum is blue-shifted by about 8 
km s~ . Perhaps, the emission at maximum light is related 
to the fact that RY Sgr is a large-amplitude radial velocity 
pulsator. There is no hint of emission in our spectra of the 
other RCB stars, even in V482 Cyg observed below maxi- 
mum light. 

By inspection, the 9850 A blend of Fe 11 and [C 1] lines 



is very similar to that of XX Cam and, hence, the forbid- 
den line is a minor contributor to the absorption feature. 
Synthesis of the 8727 A line suggests that emission overlays 
the absorption line and one might suppose that an abun- 
dance log e(C) ~ 9.2 is a lower limit but this assumes that 
the emission comes from a region exterior to the regions in 
which the absorption line is formed. This is not necessarily 
a valid assumption for RYSgr with its large pulsation and 
evidence for an internal shock. We decline to quote a car- 
bon abundance. But, a clue might be derived from RYSgr 
observations that the emission effects both [C i] lines. The 
estimations from the red absorption wing (assumed to be 
unaffected by emissions) of both these [C i] lines suggests 
the same carbon abundance of about 9.2 which agrees with 
that from C I lines. 

With the adopted model, the CN lines return a N abun- 
dance about 0.6 dex higher than from the N I lines used here 
and by Asplund et al. This difference between N abundances 
can be eliminated by lowering the adopted T e g by about 200 
K, or by a combination of a smaller T B a change accompa- 
nied by a reduction of logg. These adjustments are within 
the uncertainty limits considered reasonable by Asplund et 
al. 

5.4 UV Cas 

The 8727 A line is reasonably well fit with the abundance 
loge(C) = 8.3 (Figure 4). Two observations show no con- 
vincing evidence for the 9850 A [C i] line, indicating log e(C) 
< 8.7 (Figure 3). The two Cl lines in the 9850 A region 
suggest an abundance close to log e = 9.0, a value only 0.2 
dex less than that obtained by Asplund et al. (2000). The 
CN lines are not detectable in our spectra. UVCas joins 
XX Cam in showing a larger carbon problem for the 8727 A 
line than for the C I lines. 

5.5 VZ Sgr 

This RCB star was classified as a 'minority' star by Lambert 
& Rao (1994), i.e., the metal lines are weak. This is evident 
from Figure 1 which shows that the Fel and Fell lines are 
absent from VZ Sgr's spectrum. Unfortunately, the region 
around 8727 A has not yet been observed. Synthesis (Figure 
3) suggest an upper limit log e(C) < 8.7 for the [C i] 9850 A 
line, a limit clearly below the input abundance of 9.5. The 
two C I lines give log e(C) =9.1. Asplund et al. gave log e(C) 
= 8.8 ± 0.3. 

5.6 R CrB 

At 8727 A the synthesis corresponding to log e(C) ~ 8.4 
provides a fair fit to the observed high-resolution profile of 
2002 July 31 (Figure 4). A lower resolution (R = 60000) 
spectrum of 1995 June 17 provides the same abundance. In 
neither spectra is there a hint that emission has distorted 
the absorption profile. 

The 9850 A syntheses include the CN lines blending 
with the [C i] and Fe II feature. On the high- resolution spec- 
trum of 2002 July 31, the CN 1-0 lines are blue-shifted by 
3.6 km s _1 relative to the high-excitation lines of Cl and 



N I. Since the [C i] line and CN lines are all of low excita- 
tion, we assume that they have a common blueshift. Syn- 
thetic spectra (Figure 3) indicate an abundance loge(C) = 
9.1. Examination of our library of 9850 A spectra of RCrB 
(e.g., Rao & Lambert 1997; Rao et al. 1999) shows that the 
velocity shift between CN and the C I and N I varies during 
the star's pulsation and this shift ranges from —3.0 to +4.0 
km s _1 . Syntheses taking into the account the velocity shift 
return the same abundance from four different spectra. If 
the blue-shift is neglected, the fit to the 9850 A feature is 
less satisfactory but the derived carbon abundance is little 
affected. Asplund et al.'s atmospheric parameters are used 
in all cases. The C I lines near 9850 A are well fitted with a 
similar abundance (loge(C) = 9.2), an abundance equal to 
that determined by Asplund et al. 

The 9850 A line confirms the abundance obtained from 
C I lines, but the 8727 A line gives an abundance 0.7 dex 
below that from the 9850 A line or 1.1 dex below the input 
carbon abundance of the model. 

The N I lines in the 8727 A and 9850 A regions give the 
abundance loge(N) = 8.2, a value consistent with the result 
of 8.4 ± 0.2 given by Asplund et al. The observed CN 1-0 
lines are well matched with the abundance log e(N) = 8.2. 
This may be fortuitous agreement because the sensitivity 
of the CN line strengths to a change of the atmospheric 
parameters, especially to T e g, is high. Note that a change of 
250 K changes the required abundance by 0.5 dex. 

5.7 SU Tau 

The abundance log e(C) = 8.2 provides an excellent fit to the 
observed [Cl] line at 8727 A (Figure 4), whereas loge(C) = 
8.9 is required for an equivalent fit to the observed 9850 
A feature (Figure 3). For this fit, the CN lines are blue- 
shifted by 6 km s^ 1 from the high excitation atomic lines. 
This shift is assumed to apply to the [C i] line also, but the 
derived abundance is not critically influenced by the shift. 
The abundance from the 9850 A [C i] line is consistent with 
that derived by us from the 9850 A region's C I lines and the 
Ci lines used by Asplund et al. In contrast and consistent 
with the results for XX Cam and RCrB, the 8727 A line 
gives a markedly lower carbon abundance. 

The N abundance reported by Asplund et al. was based 
on just two N I lines. Examination of our superior spectra 
provides a more accurate equivalent width for one line. The 
other line is not present on our spectrum. We identify a 
third line. We adopt Wiese et al.'s (//-value and derive the 
N abundance loge(N) = 7.9. We suggest that this is a more 
reliable estimate than Asplund et al.'s value of 8.5. The N 
abundance estimate from the CN lines is 0.5 dex lower. The 
difference is erasable with minor adjustments to the adopted 
atmospheric parameters, such as an increase of T e ft by only 
140 K. 

5.8 V482 Cyg 

There is a hint of emission at the 9850 A [C i] line. Exam- 
ination of the three individual exposures shows that this 
'emission' occurs in only one exposure that being the one 
with the lowest signal-to-noise ratio. We, therefore, consider 
the emission to be an artefact which is removed in Figure 3. 



A fit to the absorption feature suggests an abundance upper 
iimit log e(C) < 8.9. 

In conflict with this upper limit, the C abundance of 
loge(C) = 9.5±0.3 is derived from Cl lines in our spectrum. 
Asplund et al. obtained log e(C) = 8.9 from C I lines. Ex- 
amination of our spectrum shows that C I lines across the 
spectrum are systematically stronger than reported by As- 
plund et al. This strengthening and the adoption of the same 
model as Asplund et al. necessarily results in a higher carbon 
abundance. Recall that the star was two magnitudes below 
maximum light. We suppose that the atmospheric structure 
differed from that prevailing at maximum light. The tem- 
perature gradient may have been steeper in decline than at 
maximum light. A very similar strengthening of Cl (and 
other) lines was noted by Rao et al. (1999) in spectra of 
R CrB taken early in a deep decline. 

The Nl lines give an abundance 0.8 dex higher than 
that from the CN lines. Erasure of the N I - CN discrepancy 
requires a different choice of atmospheric parameters. The 
possibilities range from raising T e g by 400 K at the adopted 
logp to lowering \ogg by about 0.8 dex at the adopted T cS . 
These are not unacceptable changes given that the atmo- 
sphere may have been perturbed. 

5.9 GU Sgr 

The blend containing the 9850 A [C i] line is contaminated 
with CN lines which are strong in this spectrum. Our synthe- 
ses suggests an upper limit log e(C) < 9.0, a value consistent 
with Asplund et al.'s value of 8.8. 

Analysis of the CN lines gives a N abundance of log e(N) 
= 8.1. Since the N I lines in the 9850 A region are blended, we 
compare the CN-based abundance with the result log e(N) 
= 8.7 ± 0.5 given by Asplund et al. The difference may be 
removed by a modest (250 K) increase of the adopted T e g 
at the adopted \ogg. 

5.10 Summary of the Carbon Abundances 

Our analyses use the appropriate MARCS model for the at- 
mospheric parameters recommended by Asplund et al. and 
C/He = 1% by number. This input abundance corresponds 
to a carbon abundance log e(C) = 9.5. In Table 2, we sum- 
marize the carbon abundances derived from the two [C i] 
lines, and the C i-based abundance given by Asplund et al. 
illustrating the carbon problem. 

Except for V482 Cyg, the carbon abundance derived 
from the two 9850 A region C I lines is in good agreement 
with that obtained by Asplund et al. 

Table 2 shows that a carbon problem extends to the 
9850 A forbidden line. This conclusion is contrary to that 
reached by Asplund et al. from analysis of equivalent widths 
of the 9850 A line. The explanation is that Asplund et al. did 
not recognize that the [C i] line was a blend. In the case of 
V482 Cyg, which was observed two magnitudes below max- 
imum light, the carbon problem has vanished for the Cl 
lines but not for the 9850 A [C i] line. This is the only case 
where the C I lines and 9850 A differ in the abundance they 
provide. 

Surprisingly, the 8727 A [C i] line offers further infor- 
mation on the carbon problem. For each of the four stars 



for which we have observed the 8727 A line, the derived 
abundance is less than that from the C I lines. In the case of 
RCrB and SUTau, the 8727 A line gives an abundance 0.7 
dex less than that derived from the 9850 A line. That the 
two forbidden lines give very different abundances is espe- 
cially puzzling if these lines are formed, as expected, in or 
close to LTE and, as might be supposed, are a product of 
the stellar photosphere, i.e., a region with the temperature 
decreasing monotonically outward. An unidentified atomic 
line superimposed on the 8727 A forbidden line would mean 
that we have overestimated the carbon abundance from that 
line. In this scenario, the carbon problem is more severe for 
both of the [C i] lines than for the C I lines. However, it is 
difficult to find a carrier for a line which is strong in RCrB 
and SU Tau but weak in other RCB stars and 7 Cyg. 

The sensitivities of the permitted and forbidden lines to 
the choice of model atmosphere are different. To illustrate 
these sensitivities, we give in Table 3 the abundances derived 
from RCrB's lines for a series of model atmospheres centred 
on Asplund et al.'s choice of (T c a , log g, £) = (6750, 0.5, 7). 
The carbon abundance from C I lines is, as expected from 
the insensitivity of the ratio of line to continuous opacity 
to physical conditions, almost independent of the choice of 
model. The abundance is also insensitive to the choice of the 
microturbulent velocity £ for a weak line but not for a strong 
line. In Table 3, the mean C abundance derived from weak 
and strong permitted lines is given. The abundance from 
the forbidden lines is insensitive to the choice of the surface 
gravity but dependent on T e ff. We note that the change ± 
250 K in T e g leads to a change in the carbon abundance 
of ± 0.2 dex. There is no dependence on £. The T e s and 
log g sensitivities are not very different for the hottest stars 
(XX Cam and RYSgr) and the coolest (GUSgr). 

To increase the carbon abundance from the 9850 A [C 1] 
line to the input abundance (loge(C) = 9.5), requires that 
the T cB of RCrB and SUTau be raised about 500 K, and 
750 K, respectively. These increases are not only outside the 
bounds considered acceptable by Asplund et al., but they 
do not remove the carbon problem for the permitted carbon 
lines. In addition, they introduce a discrepancy between the 
N abundance from the N I and CN lines. Also note that the 
higher temperatures do not eliminate the abundance differ- 
ence of 0.7 dex from the 9850 A and 8727 A forbidden lines 
of RCrB and SUTau. In short, the [Cl] lines are part of a 
now enlarged carbon problem. 



6 DISCUSSION 

Possible solutions to the carbon problem presented by the 
C 1 lines were discussed by Asplund et al. (2000 -see also 
Gustafsson & Asplund 1996). Two proposed solutions were 
suggested by Asplund et al. as worthy of further considera- 
tion: the gf- values of the C I lines are in error, or theoretical 
model atmospheres are a misrepresentation of the tempera- 
ture structure of real RCB stars. In addition to commenting 
on these solutions as they affect the forbidden lines, we ex- 
amine departures from LTE as they affect line formation. 
Finally, we discuss hand-crafted model atmospheres includ- 
ing models with a chromosphere (a temperature rise in the 
outer layers). 



Table 2. Photosphcric C abundances from Cl and [Cl] lines 
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R CrB 


6750 


0.50 
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7.0 
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V482 Cyg 
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0.50 
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6100 


0.55 


3.5 


7.9 




<8.2 


7.9 



"Asplund ct al. (2000) for R CrB stars, the expected C abundance is 9.54 when 
C/He of 1% models arc used; Luck & Lambert (1981) for 7 Cyg 
based on two C I lines near 9850A 
c scc Section 5.1 of the text 



Table 3. Carbon abundances for RCrB 
for various model atmospheres 
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9.3 
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8.4 



6.1 The gf- values 

The gf- values of the [C 1] lines are known accurately from 
theoretical calculations. Confirmation and extension of the 
carbon problem from the forbidden lines strongly suggests 
that the gf '-values of the permitted lines can not be primar- 
ily responsible for the problem. 

6.2 Non-LTE Effects 

Non-LTE effects in neutral carbon atoms in RCB star atmo- 
spheres were evaluated by Asplund & Ryde (1996). Their 
results indicate that departures from LTE are confined to 
shallow optical depths. Effects on the C I lines are slight 
because a line and the continuum are formed deep in the at- 
mosphere between levels in the carbon atom with small and 
very similar departures from LTE. Gustafsson & Asplund 
(1996) give the correction to LTE abundance from Cl lines 
as less than 0.1 dex and suggest that +0.02 dex is a typical 
value. Thus, non-LTE effects can not account for the 0.6 dex 
typical carbon problem arising from the permitted lines. De- 
parture coefficients given by Asplund & Ryde show that a 
non-LTE correction to the LTE abundance from a forbidden 
line is also very small and cannot be the source of the line's 
carbon problem. These non-LTE calculations adopt a model 
computed in LTE. 



Carbon problems may be reduced or even eliminated by 
invoking a chromosphere. A temperature rise at the top of 
the photosphere is possible even in the absence of mechani- 
cal energy deposition. This non-LTE effect sometimes called 
the Cayrel mechanism (Cayrel 1963) is discussed by Mihalas 
(1978). To assess the temperature rise for a RCB star due to 
the Cayrel mechanism, non-LTE effects on the populations 
of the carbon atom were included in the computation of the 
model atmosphere. These non-LTE effects on the carbon 
level-populations were calculated using MULTI (Carlsson 
1986), with the modelled carbon atom described in Asplund 
& Ryde (1996). Realistic background fluxes, computed with 
the MARCS model atmosphere code, were included. Typical 
departure coefficients (b=n/n(LTE), where n is the occupa- 
tion number density) of about b = 1 for the lower C I levels, 
b = 1.4 for intermediate ones and b = 0.6 for the upper 
ones were obtained. These figures were found to be char- 
acteristic of the depth interval —4 < log tr osb < —3, with 
much smaller and opposite effects at larger optical depths. 
These departures from LTE should mainly affect the opac- 
ities. C 1 is also an important electron contributor, but our 
calculations show only an underionization of less than 10% 
at shallow depths. 

We then brought these results into the MARCS pro- 
gram in order to include their effect on the continuous opac- 
ities. This procedure is not fully self-consistent - the models 



are still LTE models - but the main effect of the departures 
on the opacities are modelled in a reasonable way. They es- 
sentially diminish the C I opacities longwards of 6000 A, to 
0.6 times the LTE value, and increase them shortwards of 
that wavelength to 1.4 times the LTE value. This was ap- 
plied for the uppermost atmosphere with a decrease corre- 
sponding to our non-LTE results at greater depth. A Cayrel 
effect resulted, as expected, with a heating of the upper lay- 
ers, but only by about 1-6 K, the maximum occuring at log 
TRoss = —2. 

Thus, the resultant Cayrel effect is small for C I in RCB 
stars with a temperature increase less than 10 K at the top 
of the photosphere. This micro-chromosphere does not alle- 
viate the carbon problem! The small Cayrel effect is partly 
linked to the fact that carbon atoms are not totally domi- 
nating the scene. In fact at these relevant depths photoion- 
ization of carbon atoms only contributes typically 50% of 
the opacity. So, even drastic changes (at least reductions) to 
the carbon opacity only lead to moderate changes of fluxes 
and of the energy balance. Another circumstance is that the 
ultraviolet flux in the surface layers is so weak that it does 
not seem to matter very much and is unlike the situation 
for hot stars and the H I opacity, only a tiny fraction of the 
flux of typical RCB stars comes at energies higher than the 
ionization threshold, or even at energies able to photoionizc 
from lower excited C I levels. 

6.3 Hand-crafted Models 

The fact that the carbon problems defined by the C I and the 
[C i] lines are very similar from star-to-star across the sample 
of analysed RCB stars implies that the problems' resolution 
cannot depend very sensitively on a star's individual charac- 
teristics. Departures of the atmospheric structure from that 
predicted by a standard (MARCS) model must be very sim- 
ilar across the sample of stars. If departures are attributable 
to deposition of mechanical energy, the flux of mechanical 
energy cannot vary widely across the sample. 

That the atmospheric structure is partly the culprit is 
suggested by the observation that an RCB star in a decline 
shows a reduced carbon problem. The case of V482 Cyg was 
discussed in Section 5.8. There, the carbon problem for the 
C I lines vanished. Note that the C I lines are noticeably 
stronger in our spectrum than in spectra taken at maxi- 
mum light. The problem remains for the 9850 A forbidden 
line, but this might be the result of overlying emission, as oc- 
curred in the case of RY Sgr. (Unfortunately, the 8727 A line 
was not on the recorded spectrum.) A similar strengthening 
of C I lines leading to elimination of the carbon problem but 
with retention of the problem for the [C i] 9850 A line was 
seen by Rao et al. (1999) for RCrB in decline. 

6.3.1 The photospheric temperature gradient 

The carbon problem may be related to an incorrect rep- 
resentation by the MARCS RCB star models of the tem- 
perature gradient through the region of line formation; a 
steeper/shallower gradient produces a stronger/weaker line. 
In this connection, it may be noted that Lambert & Rao's 
(1994) abundance analysis used unblanketed model atmo- 
spheres (Schonberner 1975) and found fair agreement be- 
tween the observed and calculated equivalent widths of the 



C I lines. Note that the small carbon problem of 0.2 dex 
might be attributed to a combination of identifiable uncer- 
tainties. The primary factor responsible for the negligible 
carbon problem is that Schonberner's unblanketed models 
have a shallower temperature gradient in the line forming 
region than the MARCS (line blanketed) RCB star models 
(Gustafsson & Asplund 1996). 

Unblanketed models were calculated with the MARCS 
code. Our abundance analysis for R CrB using unblanketed 
MARCS models shows that the carbon problems presented 
by permitted and forbidden C I lines are similar for unblan- 
keted and blanketed MARCS models. The carbon abun- 
dance derived from [C i] 9850 A line using unblanketed mod- 
els again agrees that from the C I lines. The large difference 
between the carbon abundance derived from the [C i] 9850 A 
line and the [C i] 8727 A line using the blanketed models is 
unaltered by using the unblanketed models. The difference 
between the unblanketed MARCS and Schonberner's mod- 
els reflects differing temperature gradients of the models. 

Asplund et al. describe hand-crafted adjustments to the 
line blanketed models that reduce the photospheric temper- 
ature gradient and alleviate the carbon problem presented 
by the C I lines. A shallower gradient (with respect to a blan- 
keted MARCS model) was qualitatively justified as arising 
from deposition of mechanical energy whose source may be 
the sub-photospheric convection zone. Examples of temper- 
ature modified models, which alleviate the carbon problem 
from C I lines, were kindly provided by Martin Asplund (pri- 
vate communication) . Our analyses using these models show 
that the carbon problem vanishes for the 9850 A [C i] line as 
it does for the C I lines, but the 8727 A [C i] line continues 
to present a problem; the carbon abundance from 8727 A 
[C i] line is about 0.7 dex less than that from 9850 A [C i] 
line, a difference found also for unmodified line blanketed 
and unblanketed models. 

6.3.2 A chromosphere 

Addition of a temperature rise ('a chromosphere') at the top 
of the photosphere may also modify the predicted strengths 
of the absorption lines in the emergent spectrum. The chro- 
mosphere must be crafted by hand without theoretical guid- 
ance, but mindful of the discussion in the previous section. 

Experiments were made by adding a chromosphere to 
a MARCS line-blanketed model. Figure 5 shows one exper- 
iment appropriate for RCrB itself. A temperature rise is 
introduced at logTR OSS = —1.6. A flat temperature profile 
of 7000 K in the interval logTR OSS < — 1.6 and a LTE spec- 
trum synthesis reproduce the profiles of the observed [C i] 
9850 A and [C i] 8727 A absorption lines. The two lines now 
return the same C abundance which is equal to the input 
abundance, or, the carbon problem vanishes for these lines. 
The weaker observed C I lines are also reproduced with the 
input abundance. Very strong C I lines are predicted, as ex- 
pected for an LTE calculation, to have emission cores. The 
emission cores might be reduced if the non-LTE effects were 
taken into consideration. Quite obviously, the chosen chro- 
mosphere is not a unique solution. If the onset of the sharp 
temperature rise is placed as in Figure 5, a chromospheric 
temperature of 6750 to 7000 K provides acceptable solu- 
tions to the carbon problems. Similar chromospheres remove 
the carbon problems for the other RCB stars. The chromo- 



standard model, Asplund et al. 2000 
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Figure 5. Three T — tr oss structures are shown - see key at the 
top of the figure. The solid line shows a MARCS model for T e g = 
6750 and logg = 0.5. The modified T — r structure with a heated 
plateau devised by Asplund ct al. (2000) to eliminate the carbon 
problem for C I lines is shown by the dotted line. Our model with 
a chromosphere is shown by the dashed line. 

spheric temperature may scale with a star's effective tem- 
perature, but the range of sampled effective temperatures is 
small. 

Addition of a chromosphere not only changes the pre- 
dicted strengths of the C I and [C i] lines but also the pre- 
dicted strengths of other lines and, therefore, the derived 
abundances. High-excitation lines such as the N I line at 8729 
A are so weakened by chromospheric emission that the ob- 
served absorption line in XX Cam and R CrB cannot be fit 
whatever the adopted nitrogen abundance. Although careful 
crafting of the chromospheric temperature profile may alle- 
viate this difficulty, it seems likely that a non-LTE synthesis 
may be necessary. One may need to adopt spherical rather 
than a plane-parallel geometry. The chromosphere, if it ex- 
ists, is likely to be highly-structured and a far cry from the 
homogeneous plane-parallel (or spherical) layers assumed for 
the models. 

Although the artifice of a chromosphere serves to reduce 
or even eliminate the carbon problems, its presence presents 
another puzzle. The energy flux needed to produce the extra 
heating suggested here is substantial. Following Asplund et 
al. (2000), we find it to be on the order of 10 - 15 % of 
the total stellar flux. For comparison, the chromospheric- 
coronal heating of normal late-type stars is limited to about 
1 % of the total bolometric flux. A clue to the problem may 
lie in the very large widths of RCB star lines - compare the 
widths of lines (Figure 1 and 2) in 7 Cyg and the RCB stars. 
These large widths imply a very turbulent atmosphere and 
dissipation of energy. The possibility of an external radiation 
field generating the very large RCB star line widths needs 
to be explored. 

Emission lines appear when an RCB star goes into de- 
cline. The emitting gas is not necessarily to be identified 



with the hand-crafted chromosphere. Emission from a shell 
around the star would dilute the photospheric absorption 
lines. Emission in the 9850 A and 8727 A lines might weaken 
the absorption lines such that the emission is not seen but 
the weakened absorption lines return different carbon abun- 
dances. The [C 1] lines are seen in emission in spectra taken 
in decline (Rao & Lambert 1993; Rao et al. 1999). The mea- 
sured fluxes in decline are possibly lower limits to the fluxes 
at maximum light as the dust cloud responsible for the de- 
cline presumably obscures parts of the emitting shell. At 
maximum light, the cores of the strongest low excitation 
absorption lines of singly-ionized metals appear distorted by 
emission (Lambert, Rao & Giridhar 1990; Rao et al. 1999). 
It is just such lines which are prominent in the emission line 
spectrum. The emission line spectrum is of low excitation 
and does not include the C I lines. Therefore, the emitting 
shell which may contribute to the [C 1] problem is not a 
player in the Cl problem. In this context, it is relevant to 
recall the case of V482 Cyg observed below maximum light 
where the carbon problem is absent for the C I lines but not 
for the [C 1] 9850 A line. Our speculation is that the emit- 
ting shell fills in the forbidden line and the photospheric 
temperature gradient is steeper than at maximum light. 



7 CONCLUDING REMARKS 

The intriguing carbon problem presented by the RCB stars 
was discovered by Asplund et al. (2000) in an application 
of MARCS RCB star models to analysis of C I lines. Our 
contribution to the problem has been to present and analyse 
observations of the [C 1] 8727 A and 9850 A lines, also with 
the MARCS models. The 9850 A line presents a problem 
of similar magnitude (about 0.6 dex) as C I lines, but the 
carbon problem presented by the 8727 A line is more severe 
by about 0.7 dex than the original problem (about 0.6 dex). 

The fact that the carbon problems defined by the C I 
and the [C 1] lines are similar from star-to-star across the 
sample of analysed RCB stars implies that the solution can- 
not depend on a star's individual characteristics. Departures 
of the atmospheric structure from that predicted by a stan- 
dard (MARCS) model must be similar across the sample of 
stars. If departures are attributable to deposition of mechan- 
ical energy, the flux of deposited mechanical energy cannot 
vary widely across the sample. 

Uncovering the carbon problem was a surprise. Confi- 
dence in model atmospheres was shaken. Extension of the 
problem to the [C 1] lines is also discomforting. In spite of the 
fact that the RCB stars in several ways are very complex, 
they have the one property that should make the analysis of 
C 1 lines in principle simpler than for almost all other stars. 
This as the fact that the lines are due to the same element is 
the dominant opacity source and lines and opacity originate 
from levels of similar character and excitation. This provides 
a test, that is not offered for other stars. This test, however, 
fails. 

Perhaps, the problems are restricted to very peculiar 
stars such as the rare RCB stars. But this supposition de- 
serves to be tested. Luminous normal (i.e., H-rich) super- 
giants of the temperature of the RCB stars are the stars 
from which chemical compositions of external galaxies are 
now being derived. Are their atmospheric structures reliably 



simulated by our codes? Or is there a problem analagous to 
the carbon problem awaiting discovery? 
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